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Abstract 

The S.I.R. model (Susceptible, Infected, Recovered or Died) was proposed by chemistry Willam 
Kermack (1927) and the mathematician G. Mc. Kendrick (1932). the model supposes to divide 
to the individuals of a population in three categories. Susceptible to be infected, Infected and 
Recovered (immune or died by the disease). On the other hand has been a similarity in the evolution 
of epidemics of infect aerial, the computer science propagation of virus and the propagation of social 
paradigms(fashion, rumor, etc.) it calls modernly "Social Epidemics". In this work it is tried to use 
this model in different types from social networks, real or not. In order to evaluate an meta-analysis 
of results allows to investigate under wich conditions the topology of the network is excellent and 
that networks are equivalent. The result obtained can have relevance in study of propagation of 
epidemics of different types. 

PACS numbers: 05.45.-a 
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I. INTRODUCTION 



The use of mathematical models to study the evolution of epidemics is not a newness^. 
At the moment it is known enough on behavior of transmission of diseases according to 
its characteristics like knowing the limits and the possibilities these models. The S.I.R. 
model developed by Scottish chemistry William Ogilvy Kernack (1927)^ and mathematician 
A.G. Mc Kendrick (1932)^ was one on the first attempts. Unfortunately the articles fell 
in the forgetfulness until 1979 which a well-known article of Anderson and May in the 
Nature magazine^ turned east model in the departure point of the present studies. On the 
other hand, model S.I.R. generates a description in three phases of the course of an epidemic: 
starting (of slow growth), explosive and remission. Implicitly it is assumed that the contacts 
between the members of a population are purely random. This can even give relatively fit 
results for many diseases in virus transmission in networks of computers. Where it would 
be enough that it is considered to the users like susceptible individuals to the disease. This 
is thus because in general the electronic mail or the air (in the case of influenza) are opened 
enviroment that contacts to all with all with relative promiscuity then the randomness 
hypothesis can work enough^. In this work different types from real social networks obtained 
by real or fictitious study {Scale-Free SF)^ and Small World SW)^. Where a population 
of survivors was considered soon to carry out Meta Analysis of the data and thus to verify 
under which conditions the topology of the network is excellent and to investigate if the 
fictitious networks are equivalent to the real ones. 

A. Theoretical bases 

A network is a set of relations (bonds or edges) between a defined series of elements (nodes, 
vertices or actors). Formally a network is a graph defined as the triada G = (V, E, 7) where 
Vis the set of vertices, E is the set of edges and g : E V so that g{e) = {v, w}. I.E., g(.) 
assigns to each edge a pair of vertices. Nowadays to facilitate the study of the networks the 
binary matrices instead of the whole matrices are used. Then an isomorphism f : G ^ Bn 
can be defined, where -B„ is a binary and symmetrical matrix of well-known dimension nxn 
this matrix is called adjacency matrix (AM). The sociologists take like convention that in the 
rows is located the actors (exits or egos), however in the columns are located the attributes 
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or related actors (entered or alter), this convention is used in this work. As it is possible 
to be appreciated the AM contains all sensible information on the network reason why it 
allows to elaborate algorithms in MA without considering to the network. Another form to 
characterize to the network is through the well-known histogram of rows also named prestige 
of an actor who in the case of the graphs agrees with the histogram of columns or popularity 
of an actor. At the moment the connection probability is used, P{k). P{k) is the probability 
that a chosen vertex at random has k edges. According to the functional form of the tail 
of this histogram, k —* oo, the type of network can be classified in: exponential, when 
P{k) ^ e"'^'^; scale-free, when P{k) ^ k~'^~"' with 7 > 0; broad-scale when he is scale-free 
with a steep cut; and single-scale when it has a fast decay. 

II. AN AD-HOC ALGORITHM OF EPIDEMIC 

The algorithm of epidemic developed in this work supposes to know AM. This matrix 
can be obtained from investigations in a closed community, is to say does not interact with 
other communities, like as a list of electronic mail or an isolated population. Or of an opened 
community, it is to say interacts with with other communities, as it is a system of the news 
in Internet or countries connected by massive means of transport of people. To each actor 
has a structure of data with three following properties is assigned: 

1. Susceptibility or threshold (U): It is the minimum value from which the actor 
interchanges information with its pairs. 

2. Internal state: It is a state represented by a binary vector of 32 bits. 

3. Health: It can have three possible states, alive or susceptible (state S), ill (state E) 
and dead (state M). 

Of course when the actor dies automatically they eliminate the connections that connect 
with a the network. Unlike model S.I.R., in this work one assumes that an actor who 
becomes ill can cure itself but he does not stop being susceptible to return to become ill. 

On the other hand the edges of the graph are characterized by the following properties: 

1. State: It is a binary value that single taking two values "0" to indicate that there is 
not connection between actor "i" with actor "j". Or "1" to indicate that a relation 
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between both exists. 

2. Infecting: It is an assigned value of random form and uniforms in the rank C G [0, 1]. 
This it indicates the probabihty of infect between actors and it is not the same one 
for each pair of actors. 

The algorithm consists of interacting all the population of connected actors (this is defined 
as epoc) considering that each pair of related actors interacts when in a Bernoulli raffle with 
probability given by value C be favorable. Soon the condition of following proximity is 
verified. Two actors are next or neighboring if the distance of Hamming between its vectors 
of state is minor that the minimum of its thresholds. Then if this condition is verified two 
actors can interact. In the following listing the pseudocode shown the algorithm of sharing 
of information between the proximity actors Ai and A2. 

Share Algorithm 
for each actor nondead I and each actor nondead J do 
U<- mnimun(threshold(I) , threshold( J) ) ; 
if Hamining.distaiice (I , J) < U then 
for each bit do 
q<- Bernoulli raffle with probability C(I,J) 
if q>0 then J. bit <- I. bit 

else q <- Bernoulli raffle with probability C(I,J) 
if q>0 then I. bit <- J. bit 
end if 
end if 
end for 
end if 
end for 
end Share 

As it is possible to be appreciated the sharing of information is probabilistic and not 
deterministic as the case of the opinion models^. However if the actors are not verified the 
following proximity condition the operation does not take place and their vectors of state 
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do not change. This criterion is of agreed with the sharing of information codified in a data 
channel. 

The criterion of heahh adopted for an actor A is related to the value of the norm of 
Hamming of its internal state. It is to say an actor it remains in state S if the norm of 
Hamming of its internal state, H{A), is greater than the threshold than this has assigned, 
Ua- However if it happened the opposite the actor turn to the state E and if after the next 
interaction it happened that H{A) < Ua the actor turn to state M soon to disappear of the 
network assigning the value "0" to all their possible connections. Of another form the actor 
recovers his state S. 

In this model one takes like control parameter the common threshold value U to all 
the actors and the number of initial actors A^. On the other hand, initially the network 
consists of a population of K infected actors (that is to say, actors who have the state E as 
initial condition). In addition the location to infected actors is uniform in all the network 
to guarantee a random state in the population. 

III. METHODS AND RESULTS 

Four types of networks, two networks obtained in experimental works and two toy net- 
works were used, the first toy model used it was a network Small- World type (SM) with 
N=4096 actors with probability of random reconexin of 10%. The second toy model was 
a network Scale- Free (SF) with N=4096 actors with 7 = 0.1, this value is agreed with 
the experimental observations in real networks. Real social networks were obtained from 
two types of different societies in Internet. A society corresponds to a closed society as it 
happens in the lists of electronic mails (LM) and the other society corresponds to opened 
societies as it happens in news groups and the forums (DM). In the empirical networks they 
are not connected for that reason each experimental network was reduced to its connected 
maximum component. 

In all the cases one test of 512 experiment by network of 2000 times where the threshold 
varying between 10 < U < 20. The used statistical estimators were the absolute frequency of 
survival, the absolute frequency of final ill, the absolute frequency of infected initial. Except 
for the last one, all these variates with probability distribution obtained empirically from 
the tests. In Table I the probability of survival based on the threshold for the different types 
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from networks can be appreciated. 



Table I: Survival Probability 



Threshold 


SM 


SF 


LM 


DM 


10 


0.989 


0.989 





989 


0.99 


11 


0.973 


0.965 





967 


0.971 


12 


0.93 


0.883 





929 


0.907 


13 


0.827 


0.758 





891 


0.791 


14 


0.679 


0.641 





821 


0.674 


15 


0.549 


0.572 





702 


0.586 


16 


0.406 


0.448 





563 


0.485 


17 


0.298 


0.35 





246 


0.231 


18 


0.192 


0.241 





246 


0.231 


19 


0.118 


0.173 





179 


0.15 


20 


0.07 


0.096 





097 


0.07 



In most of the networks for U = 15 the survival probability it is P{S) > 50%. On the 
other hand the variance of survivors (Table II), V(S), shown that is increased until reaching 
a maximum value soon to diminish. When one studied the final networks of survivors this 
increase of the variance was due to the great dispersion of final graphs is to say: by each 
successive test the final population of survivors was totally different. 

Table II: Variance of Survival 
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18 
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1684.n3 


5590.49 


1657.92 


19 


849.35 


1194.13 


4565.42 


1301.69 


20 


574.98 


7nn.n3 


2713.59 


441.93 



Soon it was applied meta-analysis between networks of survivors for values of U=15 
by two reasons. First: probability P(S) is greater than 50% what allows to have a social 
structure not so destroyed. Second: the V(S) is great this allows to guarantee a great variety 
of final societies. 

In the following tables it is the comparative results in each particular situation. DM and 
LM correspone to the analysis group of which compares them with the control group given 
by networks SM and SF. The dead and alive numerical values correspond to the mean value 
of each test. The variable Q represents the sum in a column and the variable R the sum in 
a row. 

Table III Netowrk Meta-analysis of LM vs SM 





LM 


SM 


R 


Dead 


60 


2398 


2558 


Alive 


140 


1698 


1838 


Q 


200 


4096 





Table IV Netowrk Meta-analysis of DM vs SM 





DM 


SM 


R 


Dead 


83 


2398 


2481 


Alive 


117 


1698 


1815 


Q 


200 


4096 





7 



Table V Netowrk Meta-analysis of LM vs SF 





LM 


SF 


R 


Dead 


60 


2465 


2525 


Alive 


140 


1631 


1771 


Q 


200 


4096 





Table III Netowrk Meta-analysis of DM vs SF 





DM 


SF 


R 


Dead 


83 


2465 


2548 


Alive 


117 


1631 


1748 


Q 


200 


4096 





IV. CONCLUSIONS AND RESULT ANALYSIS 

In a prospective analysis of data where the probability of death in one given network X 

Dead 



is: 



px = P{M/X) 



Q 



X 



(1) 



After from previous tables psM = 0.58; plm = 0.30; poM = 0.41 and psF = 0.60, then if it is 
defined as null hypothesis, Hq, to not exist difference in probability between two networks. 
I.E. if Hq is accepted soon is fulfilled a 95% of certainty that pi = p2 or DR = Pi — P2 = 
and both networks are equivalents. DR is known as difference of proportions whose variance 
comes given by: 

(2) 



V{DR) = Pl^^^ + P'^'-P'^ 



Qi Q2 

here it is assumed that the random varibles have statistical independence. Soon from the 
interval of following confidence: 



DR±1M^/V{DR) 



(3) 



If DR = belongs at this interval it accepts in probability Hq other wise it is rejected and 
it posible to conclude: 

1. When it was compared networks LM and DM were / = 0.02 ± 0.09. In this cituacin 
it is clear that Hq cannot be rejected and both networks are equivalent. 
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2. When it was compared networks SF and SM were I = 0.02 ± 0.021. Marginally it 
is possible to be accepted that both networks are equivalent, but to affection of the 
analyses they can be taken like equivalent. 

3. When it was compared LM with SM were / = 0.28 ± 0.06. That clearly shown that 
there is to reject the null hypothesis. 

As a conclusion it is possible to be appreciated that there is no equivalence between 
the theoretical networks and the obtained ones by experimental work when the proposed 
algorithm is applied to them. It is to say in this prospective analysis the exposed actors to 
a network type DM or LM they have more probability of life then the behavior is diffent on 
surviving population of the found one in toy networks. This fact brings like consequence the 
nonreliability of the simulations on propagation of aerial epidemics of virus in theoretical, 
letting a door opened to find networks theoretical more realists. 



* Electronic address: lhcaste@ifir.edu. axl 
^ Electronic address: lili@ungs.edu.ar 
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